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The concept of QCD sum rules is extended to bound states composed of 
particles with finite mass such as scalar quarks or strange quarks. It turns 
■ out that mass corrections become important in this context. The number 

of relevant corrections is analyzed in a systematic discussion of the IR- and 
UV-divergencies, leading in general to a finite number of corrections. The 
results are demonstrated for a system of two massless quarks and two heavy 
scalar quarks. 
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In the last fifteen years the concept of QCD sum rules invented by Shifman, Vain- 
shtein and Zakharov |]J has been very well developed. The intention was to include 
non-perturbative effects in QCD calculations, which naturally become important 
when dealing with hadron physics. The concept of vacuum condensates is best 
illustrated in the case of chiral symmetry, which is spontaneously broken for hadrons 
Ph! as one can see by comparing the masses of different mesons with identical quantum 

numbers (i.e. p and ai meson). A direct consequence of spontaneous breakdown of 
any symmetry is the appearance of non-vanishing vacuum condensates. The idea 
of QCD sum rules is to parametrize the non-trivial QCD vacuum with vacuum 
condensates, including the chiral condensate < qq >, and in this way to realize a 
parametrization of non-perturbative effects in QCD. 

The condensates appear as matrix elements in the operator-product-expansion 
||] (OPE) of a physical function (i.e. the polarization function): 



n(g 2 ) ~ E < > / d * x E A X ) eiqx 

3=1 



This expansion separates perturbative (Wilson coemcents Ej(x)) and non-pertur- 
bative (condensates < 0|5j(0)|0 >) parts, being thus completely calculable with 
perturbative Feynman rules. It can be related to the hadronic spectrum through a 
dispersion relation with one subtraction C 



^-OPE{q ) = — ds — 7 ^~ + C (2) 

7T Jo s{s — q 2 ) 
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where the imaginary part of the function Hho4 is connected to the hadronic spec- 
trum. In that way it is possible to get hadronic quantities like the mass of a bound 
state as a function of a few vacuum condensates, which are fitted to a large number 
of hadronic data. 

The kind of condensates being relevant in the OPE depends strongly on the 
flavour composition and the quantum numbers of the bound state considered. In 
the case of light quarks (up or down) there is strong evidence for the presence of 
quark-antiquark pairs in the vacuum, so that apart from the gluon-condensates the 
quark-condensates have to be taken into account. If the particles are heavy quarks 
(charm or heavier) the value of the corresponding condensates becomes very small 
and may be neglected, so that the calculation may lead to a reasonable result by 
including the gluon-condensate only. The strange quark, however, poses a problem 
because it is neither light nor heavy. Therefore its condensate can not be neglected 
a priori, and in order to take it into account, it is necessary to treat massive quarks 
in QCD sum rules — a problem unsolved up to now. We derive several results 
necessary to approach this goal and elucidate the problems occuring. 

This article is organized as follows. We shall illustrate the problem of finite 
mass sum rules with a specific example (introduced in section 2), a detailed study 
of which will be published elsewhere. In the third section we give an expression 
for the matrix elements of normal-ordered, nonlocal field products — which occurs 
in the expansion of the correlator under consideration — , demonstrating that the 
finite masses cause corrections to the commonly used expressions. Afterwards some 
of the Wilson coefficients introduced above will be discussed with special attention 
to the appearing infrared and ultraviolett divergencies and to the relevance for the 
sum rule. 



2 The system under consideration 

For studying the influence of mass terms in QCD sum rules calculations we will 
consider a four-particle bound state composed of two light fermionic quarks and 
two heavy scalar quarks. The last ones are hypothetical particles introduced in 
QCD phenomenologically. Scalar quarks are predicted by the GUT-theory super- 
symmetry || (they are called squarks), such that we are not merely dealing with 
a toy model. The physical problem motivating this model is wether or not the 2 
quark - 2 squark bound state may be energetically more favourable than a 2 squark 
state, which is predicted to have a very high mass |4[]. The lowest lying 2 quark - 2 
squark bound state can be calculated — using QCD sum rules — in dependence of 
different new condensates and of the squark mass. These parameters are unknown, 
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but if one could find a reasonable choice of them, such that the 2 quark - 2 squark 
bound state mass becomes small, this would be most interesting for supersymmetry 
phenomenology. The results of these calculations will be published elswhere. 

In the following calculation the quark masses are neglected and the scalar quarks 
are treated as massive particles. In order to calculate the mass of the lowest lying 
bound state of this system, the OPE of the polarization-function corresponding to 
the diagram in Figure [I] has to be considered (see eq. (|TJ)), which is given by the 
two-point-function 

m 2 ) = ~ j d\x - y) e ik ^ 9lw < 0\T{J»(x), r(y)}\0 > (3) 

where T denotes the time ordered product. The current J^(x) contains an incoming 
and outgoing fermionic and scalar quark: 

J M (x) = g : ip(x)j' x ip(x) ^(x)(j>(x) . (4) 

The present sum rule is calculated in lowest order of QCD perturbation theory. 
Thus the exchange of gluons and the coupling to the gluon condensate are neglected. 

Usually QCD sum rule calculations are carried out in configuration space. In 
the case of massless particles this choice simplifies the calculations: the propagators 
are as simple as in momentum space, but the number of integrations is reduced to 
one, while in momentum space the number of integrations is equal to the number 
of loops in the corresponding Feynman diagram. However, in the case of particles 
with non- vanishing masses the propagators become very cumbersome, so that in 
most cases calculations in momentum space will be more convenient. 

The next step is to find the explicit form of the OPE for the polarization func- 
tion (3). The time-ordered product has to be expanded into normal-ordered prod- 
ucts with attention to the non- vanishing vacuum condensates. In this way one gets 
sixteen terms, each of them corresponding to one diagram: first there is the per- 
turbative diagram (see figure P, which corresponds to the fully contracted term 
without any condensates 

- < | Tr{S(x - y)^S(y - x)Y} Tr{A(y - x)A(x - y)} \ > (5) 

where S(x — y) is the Feynman-propagator of the quarks and A(x — y) is the scalar 
quark propagator. The traces are taken over color- and Dirac-indices in the case of 
fermionic propagators and in color space only in the case of scalar propagators. 

There are eight terms containing one non-contracted pair of fermionic field 
operators of the form 

+ i < | (Y) lk S kl (y - x) ( 7 % : ^(y^x) : Tr{A(y - x)A(x -y)}\0> (6) 
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Figure 1: Perturbative Feynman diagram with two scalar (normal lines) and two 
fermion propagators (heavy lines) 

which all vanish independently of the scalar part of the term: the condensate of 
massless fermions is proportional to 5ij, such that a trace over an odd number of 
7-matrices appears in the fermionic part of the expression. 

The term without any contraction is vanishing as well, because it contains no 
propagator, so that there is no momentum flow through the diagram. 

Six diagrams are left with one or two non-contracted pairs of scalar field oper- 
ators, for example the contribution from figure |2| 

- i < I Tr{S(x-y)^S(y-x)Y} A{x-y) : 4>{x)<j>{y) :| > . (7) 

When expanding expressions like that given in eq. @ one quite generally encounters 
nonlocal normal ordered product of scalar or fermionic field operators. In these 
products the field operators are taken at different points in Minkowski space, while 
the vacuum condensates of scalar or fermionic fields are defined at the same points. 
Nevertheless, the nonlocal normal ordered products commonly are identified with 
the vacuum condensates. We claim that this procedure is correct only in zeroth 
order in the mass of the particles under consideration. 



3 Mass corrections to the condensates 
3.1 Condensates of scalars 

In the previous section we have shown, that vacuum expectation values of normal 
ordered products of scalar fields at different space-time points of the type 

<O|:0 a (a;)^(y):|O> (8) 
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occur in the operator product expansion. Here the greek letters denote the color 
space indices important in the following. With the object of establishing a relation 
between the expression (|8|) and the vacuum condensates defined by 

< # > := < °l : 0*(x)M x ) : 1° > ( 9 ) 

a/3 

we will have to expand one of the fields in eq. @ with respect to the space-time 
point y belonging to the second field. The corresponding Taylor series is 

OO 

= E - vT 1 ■■•(*- vT n {4> d, n )(y) . (10) 

To preserve gauge invariance the ordinary derivatives must be replaced by the 
covariant ones. 

OO -| _ _ 

= E-( x -vT 1 ---( x -yr n ^D^---D^)(y) (11) 

n=0 U - 

The replacement is allowed if there exists a gauge for which the above two expres- 
sions are equivalent. The coordinate gauge x^A ^ = has exactly this property (see 
for example ||). 

Using this expansion we will derive an expression of the full series describing 
the nonlocal field product in eq. (§): 

< 0| : ^ a {x)<j)p{y) : |0 > 
= <O|:0 o (y)^(y):|O> 

+ <o\-.Ud, 1 ) (y)My)-\a> ^ 



+ \ < 0| : (V 5 M1 5 M2 ) Q Q/)^) : |0 > 
+ -•• 

+ < 0| : U 5 M1 ■ ■ • 5„ m ) (2/)^(y) : |0 > e Ml " " ■ £ Mm 

+ ••• (12) 

where £ := x — y. 

The terms in eq. ([12]) can contribute to the series only if the corresponding 
operators do not have any uncontracted Lorentz index, for the vacuum state is 
a scalar state. Consequently only the terms with an even number of covariant 
derivatives are of interest (the remaining terms will be renumbered using n := m/2). 
All color and Lorentz indices must be contracted in this terms, so that the structure 
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of the vacuum expectation value in color and Lorentz space is fixed a priori: 

< o| : ~$ a (y)My) : |0 > = c <W 
< o| : (^b^b^(y)My) : |o > = d 

1 M2 9fJ,3 + 9fJ,l V3 9 ^2^4. + 9 VI HA ^2 ) <W 

< 0|: ( 0D W • ■ • B M2n ) (y)My) : l° > = C n (g^ 2 g^ 4 ■ ■ • ^ 2ft _ 1M2n + perm.) 5 Q/3 

(2n-i)H terms 

(13) 

The coefficients C n are to be calculated by contracting the equations ( Jl3| ) with the 
metric tensor g^ 1 ^ 2 , using the Klein-Gordon field equation (□ + m 2 )<f)(y) = with 
□ := g^D^Dj;, and summing over all remaining indices. By this procedure one 
gets the condensates (|9|) multiplied with the mass in the corresponding dimension 
on the left hand side and the coefficient of interest with a numerical factor on the 
right hand side. For instance: 

<#> = J2<°\-0*(y)My)--\Q> 

a/3 

= C ^25 a/3 = N C C 

a/3 

=► Co = ^-<#> (14) 
- m 2 < # > = £<0|: Ug^ B M1 B M2 ) : |0 > 

= CxY^Sap^g^ = W C X 

a/3 
2 

777 

=► d = -^<00> (15) 

a/3 V /a 

~~ C*2 fia/39^ 9^ {d 1*11^9 f>3IH ~l~ 9fJ,ifJ,-i9fi2lJ.i ~l~ 9h\Ih9hwvs) 

= N C C 2 (16 + £ + <^ \g r) = 24iV c C 2 



M3^M2 1 J M4^A*2 

,4 



ft = m <#> (16) 



where N r is the number of colors. 
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The general form of C n> 2 is 



2n 

m 



Cn = (-l) n ^-<#> , (17) 

with 

G n ■= g^ 2 ■ ■ ■ g^-^ (g^ 2 ■ ■ ■ g^ 2n _ W2n + permutations) . (18) 

For each contraction of two covariant derivatives one gets a factor — m 2 coming 
from the Klein-Gordon field equation, all terms contain the same factor N c from 
the summation in color space, and the factor G n is appearing by construction. One 
can easily verify that the relation G n = 2 n (n + 1)! holds for all n. 

In this way the coefficients C n of eq. ([T7]) are known explicitly and thus we have 
found an expression for all the vacuum expectation values of eq. ([13]). Inserting this 
into the expansion ([12]) we get: 

_ 00 1 m 2n _ 

< 0| : : |0 > = U-i r — 2nN ^ + i)l < * > Sa0 

S v ' 

(2n-i)M terms 

(19) 

The (2n — 1)!! terms contract the 2n dimensional tensor £ m • • • £ M2n yielding (2n — 
l)!!(£ 2 ) n , so that we are left with the full expression of the nonlocal field product 
in eq. (Bp to all orders in the mass: 

1 - 00 1 /m 2 f 2 \ n 

< 0| : te^M : 1° > = W Jap < 00 > E(-l) W w , (w + 1) , [~T- ) > (2°) 



The sum can be performed giving: 

< 0| : 4 a (x)Mv) ■ 1° > = TT<W < 00 > -iTT J i MZD (21) 

N c m\t\ 

This result is not so surprising, because there is a connection between the scalar 
propagator and the Bessel functions. The left hand side of eq. ( pT]) has the structure 
of a scalar propagator as well. In this sense it may be meaningful to call the 
vacuum expectation value of a nonlocal normal ordered product of two fields a 
non-perturbative propagator. 

In the case of vanishing mass the series reduces to the commonly used relation 
< 0| : a (x)Mv) ■ 10 > = ^-Sofi < 00 > (22) 
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For norivanishirig mass the higher-order terms of the series have to be taken into 
account. The question arises wether the series is convergent and which order still is 
important in a specific calculation? The answer depends strongly on the value m 2 £ 2 
where £ 2 is the distance between the two vertices at x and y, which is inversely 
proportional to the mass M 2 of the incoming state to be evaluated by QCD sum 
rules. So the ratio m/M of the composite particle mass and the bound state mass 
determine the importance of the higher-order terms. For instance m/M = 1 leads 
to a 10% correction in first order, while m/M = 2 blows up the first order correction 
to 50% and the second order term makes another 10% contribution. 

To conclude: the higher-order corrections may be neglected for systems where 
a heavy bound state is constructed from light particles. However, the higher-order 
corrections become dominant if the bound state is lighter than its constituent 
particles. For instance this could happen in the supersymmetric system introduced 



in the second section. In the last case the exact form (|20|) or (pi ) has to be used 



3.2 Condensates of fermions 

A treatment analogous to eq. (f|) is possible also for fermion fields instead of the 
scalar fields used in the previous section. This is especially important if the in- 
terest lies in the treatment of finite masses in QCD sum rules emerging in mass 
calculations of strange baryons ||. 

The following nonlocal normal-ordered product of fermion field operators is to 
be expanded 

<0|:^ fa (z)^- /9 (y):|0> . (23) 

Here the fields carry a Dirac index (latin letter) in addition to the color index 
(greek letter). The fermion condensate is defined by 

< ^ > := £ < 0| : tfMMx) ■ 10 > (24) 

ija/3 

and the corresponding expansion of one fermion field preserving gauge invariance 
reads (see eq. (|TT|)) 

ip ia {y) = x r • • • (y - x T n ■ ■ ■ w ( 25 ) 

This series can be evaluated in close analogy to the case of the scalar field treated 



in sect. 3.1. There is one important difference: In the scalar case only the terms in 
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the expansion of the field with an even number of derivatives give a contribution 
to the series (see discussion after eq. (|12"D ), while in the case of fermion fields the 
terms with an odd number of derivatives are non- vanishing, for the Dirac equation 
{i Jp — m)ip = is of first order in the derivatives. 



The Lorentz structure of the even terms is the one shown in eq. (13) with an 
additional factor 5ij. Each odd term has the same structure as the corresponding 
even term with one more derivative, except for the replacement of one by 7^. 
It turns out that the coefficients C n {n 6 0,1,2,.. .) corresponding to eq. fll7|) are 

c » - { -M < * > (26) 

where the even coefficients G2 n coincide with those for scalar fields while the odd 
coefficients are given by G2n-i = Gin- This means that the coefficient of a term 
with an odd number of derivatives is equal to the coefficient of the following term 
with an even number of derivatives (except a factor im), which are determined by 
G 2n = 2 n (n+l)! like in the case of scalar fields. After some algebra the whole series 



corresponding to eq. ( |20| ) reads: 



^ > (-i)"^feji o^r (27) 



n-l 



This result can be verified by acting with (+i $ + m) on the series ( p0[ ) for scalar 
fields 

< |: A a (x)i> jl3 (y) :| > — (i $ + m) < |: a (x)Mv) -\ > (28) 
with the replacement 

m < 00 > — > < V> V > (29) 
necessary on dimensional grounds. So that the above series can be rewritten: 

< 0| : ^ ia (x)^P jf3 (y) : |0 > 

- ^ < ^ > <&p (-l) n ^ (^) n (30) 



4/Vc 

2iV c 



Again the commonly used formula can be recovered as lowest order term in the 
mass expansion looking at eq. (|27|) 



< 0| : lU*)lMv) ■■ |0 > = < ^ > , (31) 



4iV c 
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which remains correct for vanishing mass. In the case of non-vanishing mass the 
higher-order contributions of eq. ( |30"D should be taken into account. Belyaev and 
Ioffe tried to explain the mass splitting of strange and nonstrange baryons using 
QCD sum rules 0. It would be interesting to calculate the importance of these 
mass corrections to their results. 



4 IR-divergencies and higher-order mass terms 



The relevant equation in QCD sum rules calculations is a dispersion relation with 
one subtraction (see eq. (|^)), which gives a relation between the polarization func- 
tion on the quark level — expanded with the use of the OPE (eq. ([!])) — and the 
hadronic spectrum. The sum rule is obtained by taking the Borel transformation 
of eq. (0), which has the property of isolating the lowest state of the hadronic spec- 
trum (the calculation of moments is an alternative method of isolating the lowest 
state 0). The further the second state in the spectrum is separated from the lowest 
state, the better the isolation works. A typical formula giving an explicit value for 
the mass of the lowest state is 



Jq S ° ds s ImU OPE (s) e~ 



St 



f°°ds ImU p E (s)e- st 
where t is the Borel transformation parameter. 



M 2 (32) 



The main point in this formula is the sole appearance of the imaginary part 
of Hope, so that one needs only to calculate the terms being relevant for the 
imaginary part. The typical structure of a dimensional- regularized Feynman graph 
G(q 2 ) (e = 4 - d) is: * 

G(q 2 ) ~ | ~ + const. + In [~^-J } P(l) (33) 

where P(q) is a purely polynomial expression and [i is the renormalization param- 
eter. Remembering the analytic continuation to q 2 > the whole imaginary part 
of such a graph is hidden in the logarithm: 



47T/i 2 



q 2 



+ w6(g 2 ) (34) 



In the previous section a complete expression for nonlocal field products to all 
orders in the mass was derived and the notion of a nonperturbative propagator was 
introduced. In the following we show that in general it is justified to restrict the 
calculation of a Feynman diagram containing a condensate to a finite number of 
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higher-order corrections in the mass. To illustrate the argument we treat the case of 
one condensate of scalar fields, which represents a large class of graphs of the same 
structure. The corresponding graph is shown in figure 0. The broken line represents 




Figure 2: Feynman diagram with a scalar propagator, a condensate of scalar fields 
and an arbitrary additional part 

the condensate, the full line is a scalar propagator and the hatched area represents 
an arbitrary additional structure. This could be a massless fermion loop, for exam- 
ple, which would occur in the sum rule of the 2 quark - 2 squark system introduced 
in the second section and which corresponds to eq. (0). The corresponding algebraic 
expression G{k 2 ) including the mass corrections to all orders reads 



G(k 2 ) ~ <00> - g ) 



(2tt)< 



where the Fourier transformation of eq. fl20|) has been used and the symbol D p = 
Vp)-(^-,— Vp) was introduced. The function / represents the hatched area 



V dpo ' P J \ dpo 

in figure |2| and A(p — q) is the scalar propagator. 

To evaluate the integral I p it is useful to reflect once more on the significance 
of the notion of the nonperturbative propagator A given by the series (p0[). The 
nonperturbative propagator satisfies the free Klein-Gordon equation so that the 
relation 

/ d 4 pg(p 2 ,p,. . .) A NP (p) = J d 4 pg(m 2 ,p,. . .) A NP (p) (36) 

holds for an arbitrary function g(p 2 ,p, . . .) in lowest order of perturbation the- 
ory. Note that the calculations already before were restricted to lowest order of 
perturbation theory, so that this is no additional restriction. By applying the 



replacement rule fl36|) to the scalar propagator in eq. (|35|) , it can be reduced to 
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A(p — q) = [q 2 — 2pq + ie] . The n-th order differentiation of the distribution 
5(p) is carried out by n differentiations of the simplified propagator and after some 
algebra eq. fl35|) reads: 

G( P, „ < * > / _g, , ((i _ | _<ML_ (37) 

A similar argument can be repeated for any other function g instead of g(p) = 
A(p — q). In particular a fermion propagator would lead to a similar result. 

The infrared and ultraviolet behaviour of eq. Q57| ) is determined by the lowest 
and highest power (mi and rrih) of q in the function / ((k — q) 2 ). So G(k 2 ) is UV- 
divergent if n < mfc 2 +2 is fulfilled, while it is IR-divergent if n > m 2 +2 is fulfilled. 
One important conclusion can be drawn: only a finite number of terms containing 
UV-divergencies exists. There is a problem region m ' 2 +1 < n < m ft 2 +1 with both 
UV- and IR-divergencies occuring, and finally an infinite number of terms with 
possible IR-divergencies exists. 

To illustrate the above rules, the 2 quark - 2 squark system introduced in the 
second section is treated explicitly. In this case the function / ((/c — q) 2 ) in eq. fl3"7| ) 
becomes the regularized massless fermion loop diagram, so that eq. (j3~?D now reads: 



G{k 2 




(3f 



where all proportionality factors have been reintroduced. The constant term in 
eq. ( |3"3"D is omitted, because it does not contribute to the imaginary part of G{k 2 ). 
The above rule to analyze the divergencies occuring in this expression leads to the 
following expectation: f((k — q) 2 ) ~ (fc — g) 2 has the lowest and highest power 
mi = and m^ = 2, so that UV-divergencies occur for n < 2, IR-divergencies occur 
for n > 1, and double divergencies occur for n G {1,2}. This is exactly the result of 
the straightforward regularization of the last integral in the expression (p8|), which 
can be verified by counting the number of divergent T-functions in the limit e — > 0: 



^ 2 ) - ^<#>E o 4^(- 2 n* 




r ( 3 "!) r(-2 + | + i/ 



r(n + i)r(4-n 



(39) 



The double divergent terms for n = 1,2 lead to a divergent imaginary part, for 
the integral I in eq. ( j38|) is IR- and UV-divergent at the same time, which leads to 
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meaningless results for the sum rule. A more detailed examination of this integral is 
unavoidable. Decomposing / in three parts, in which the integrand is proportional 
to k 2 , k ■ q and q 2 respectively, one can easily verify that the only part having 
IR- and UV- divergencies is the fc 2 -term. It turns out that the IR-divergency in 
this part can be regularized by introducing a small IR-cutoff parameter a 2 into the 
scalar propagator (or equivalently using a lower integration limit), reading: 

. f A , f(k-q) 2 \ k 2 , An . 

42 = v^H^v (gw+^+i (40) 

This integral has no IR-divergent part any more, so that it is possible to regularize 
the remaining UV- divergency in the dimensional regularization scheme. The result 
to order O(e) (where e = 4 — d) depends on the cutoff parameter a 2 : 



n 2 k 2 f 4 2/i e 2l ( -k 2 

4 2 = - fn u 7 + + const. + -7 In I - 



(2tt) 4 \e 2 e y " 2 ' \IQti^ 2 ) 

+ (HMtS-Hi^)) (41) 

The essential point of this result is the exact cancellation of all terms to the or- 
der 0(1) or lower which depend on the cutoff parameter a 2 . The imaginary part 
relevant for the sum rule in this way remains cutoff parameter independent. The 
whole imaginary part of eq. (p8|) becomes finite: 



Im{G{k )} = 3(2*)° < 




k 2 

This confirms the expectation that the IR divergencies are unphysical 



0(k 2 ) (42) 



Independently of the choice of the function / in eq. fl37|) it must be possible to get 
rid of all IR-divergent terms. This becomes clear when we look back to the very be- 
ginning. The OPE is the basic formula of QCD sum rules and its idea is to separate 
the long- and short- distance effects of nonperturbative QCD. The short-distance 
effects are completely embodied in the Wilson coefficients just calculated. So any 
remaining IR-divergency occuring while calculating Wilson coefficients proves an 
admixture of long-distance effects, which indicates that the OPE did not really 
separate the two scales of QCD. Novikov et al. argued || that for the consistent 
use of OPE one has to divide the integration domain (while calculating the Wilson 
coefficients) into two parts by introducing the normalization point. In this way one 
can get rid of all IR- divergencies, for only the high momentum part remains rele- 
vant for the Wilson coefficients, like it was demonstrated in the above example. The 
remaining question is wether or not the resulting coefficients depend on the nor- 
malization point. A meaningful OPE or a meaningful sum rule can be guaranteed 
only if the result is not or weakly dependent on the normalization point. 
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Supposing that the Wilson coefficients are independent of the normalization 
point separating the scales, one has the possibility to restrict the calculation of 
the coefficients to the terms containing an UV-divergent part. And the above rule 
n < mh + 2 shows that this is a finite number of terms. A confrontation with IR- 
divergent terms is inevitable only if we calculate the double divergent terms m ' 2 +1 < 
n < mfe 2 +1 where the UV and IR divergent parts have to be separated carefully 
(for example by introducing an IR-cutoff, like it was done in the last example). 

Summarizing, there is only a finite number of mass corrections in eq. ( |20| ) to be 
included in the calculation of the Wilson coefficients, for only a finite number of 
terms is UV-divergent and the IR-divergent terms cannot contribute to the Wilson 
coefficient if — and only if — the OPE is meaningful. But this last question leads 
us to the main problem of enlarging the concept of QCD sum rules to finite mass 
particles. A particle of several lOOMeV like the strange quark introduces a new 
scale in the theory so that the separation of long- and short-distance effects in an 
OPE becomes much more difficult than in the massless case and remains a delicate 
task. 



5 Conclusions 



The application of QCD sum rules to bound states composed of particles with 
finite mass leads to a correction of the vacuum expectation value of the nonlocal 
normal-ordered field product, which does not remain simply a vacuum condensate 
(eq. ( p^) for scalar fields and eq. ( p7|) for fermions) like it does in the massless case. 
We derived higher-order corrections and the results are given in eq. fl20|) for scalars 
and eq. (|26|) for fermions to all orders in the mass . The importance of higher-order 
terms grows in the same measure as the mass ratio of the constituent particles 
and the bound state mass. These investigations are of great importance for the 
application of QCD sum rules to mesons and baryons with strangeness. 

Analysing the structure of the Wilson coefficients corresponding to the higher 
order mass corrections, we found that in general a restriction to a finite number of 
corrections leads not to any neglection of mass effects, for only a finite number of 
corrections is UV divergent. On the other hand, the number of IR divergent terms 
is unknown and unlimited. The problem of IR divergencies may be circumvented by 
introduction of an IR-cutoff provided that a seperation of short and long distance 
effects in the operator product expansion is assured to a good degree of accuracy 
and provided that the final expressions does not depend strongly on the cutoff. 
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